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Abstract 

The Next-to-Minimal Supersymmetric Extension of the Standard Model (NMSSM) with a Higgs 
sector containing five neutral and two charged Higgs bosons allows for a rich phenomenology. 
In addition, the plethora of parameters provides many sources of CP violation. In contrast to 
the Minimal Supersymmetric Extension, CP violation in the Higgs sector is already possible 
at tree-level. For a reliable understanding and interpretation of the experimental results of the 
Higgs boson search, and for a proper distinction of Higgs sectors provided by the Standard 
Model or possible extensions, the Higgs boson masses have to be known as precisely as possible 
including higher-order corrections. In this paper we calculate the one-loop corrections to the 
neutral Higgs boson masses in the complex NMSSM in a Feynman diagrammatic approach 
adopting a mixed renormalization scheme based on on-shell and DR conditions. We study 
various scenarios where we allow for tree-level CP-violating phases in the Higgs sector and 
where we also study radiatively induced CP violation due to a non-vanishing phase of the 
trilinear coupling A t in the stop sector. The effects on the Higgs boson phenomenology are 
found to be significant. We furthermore estimate the theoretical error due to unknown higher- 
order corrections by both varying the renormalization scheme of the top and bottom quark 
masses and by adopting different renormalization scales. The residual theoretical error can be 
estimated to about 10%. 
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1 Introduction 



The search for the Higgs boson and ultimately the understanding of the mechanism behind the 
creation of particle masses represents one of the major goals of the Large Hadron Collider (LHC). 
Recently, the experimental collaborations ATLAS and CMS have updated their results on the 
search for the Higgs boson. Both experiments observe an excess of events in the low Higgs mass 
range, compatible with a Standard Model (SM) Higgs boson mass hypothesis close to 124 GeV 
at 3.1cr local significance as reported by CMS [l] and close to 126 GeV at 3.5a local significance 
in the ATLAS experiment (2j. This is still too far away from the 5a required to claim discovery 
and necessitates the accumulation of further data in the ongoing experiment. The slight excess of 
events in the 77 final state signature as compared to the Standard Model expectation may hint to 
the existence of new physics. One of the most popular extensions of the SM are supersymmetric 
models (SUSY) [3j. While the Higgs sector of the Minimal Supersymmetric Extension (MSSM) [2] 
consists of two complex Higgs doublets, which lead to five physical Higgs states after electroweak 
symmetry breaking (EWSB), the Next-to-Minimal Supersymmetric Model (NMSSM) [5] extends 
the Higgs sector by an additional singlet superfield S. This entails 7 Higgs bosons after EWSB, 
which in the limit of the real NMSSM can be divided into three neutral purely CP-even, two neutral 
purely CP-odd and two charged Higgs bosons, and in total leads to five neutralinos. Although 
more complicated than the MSSM, the NMSSM has several attractive features. Thus it allows 
for the dynamical solution of the fi problem [6] through the scalar component of the singlet field 
acquiring a non-vanishing vacuum expectation value. Furthermore, the tree-level mass value of 
the lightest Higgs boson is increased by new contributions to the quartic coupling so that the 
radiative corrections necessary to shift the Higgs mass to ~ 125 GeV are less important than in the 
MSSM allowing for lighter stop masse^jand less finetuning (t]-[9] . The enlarged Higgs and neutralino 
sectors, finally, lead to a richer phenomenology both in collider and dark matter (DM) experiments. 
The latter is due to the possibility of a singlino-like lightest neutralino, the former due to heavier 
Higgs bosons decaying into lighter ones at sizeable rates or due to possibly enhanced or suppressed 
branching ratios in LHC standard search channels such as 77 or vector boson final states [9j[l0], to 
cite only a few of the possible modifications compared to SM or MSSM phenomenology. 

The enlarged parameter set in supersymmetric theories allows for further sources of CP violation 
as compared to the SM, where the only source of CP violation occurs in the CKM matrix. Hence, 
the soft SUSY breaking couplings and gaugino masses as well as the Higgsino mixing parameter /i 
can be complex. While in the MSSM CP violation in the Higgs sector is not possible at tree-level 
due to the minimality conditions of the Higgs potential, it can be radiatively induced through non- 
vanishing CP phases 11-13 . Consequently, the CP-even and CP-odd Higgs bosons mix so that the 
physical states have no definite CP quantum number any more leading to substantial modifications 



in Higgs boson phenomenology 14 . The Higgs couplings to the SM gauge bosons and fermions, 
their SUSY partners and the Higgs self-couplings can be considerably modified compared to the 
CP-conserving case inducing significant changes in the Higgs boson production rates and decay 
modes. This could allow for Higgs bosons with masses below present exclusion bounds from LEP 
and possibly Tevatron and LHC as they might have escaped detection due to suppressed couplings 
involved in the various standard Higgs search channels, which would then necessitate new search 



lr rhe bulk of the radiative corrections stems from the (s)top loops. 
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strategies [l5| . 

In the NMSSM CP violation in the Higgs sector is possible at tree-level. Though spontaneous 



tree-level CP violation in the ^-invariant NMSSM is impossible due to vacuum stability 16 



explicit CP violation can be realized already at tree-level in contrast to the MSSM. In principle, 
there can be six complex phases parametrizing the CP violation in the Higgs sector, two relative 
phases between the vacuum expectation values (VEVs) of the Higgs doublet and singlet fields 
and four phases for the complex parameters A, k,A\,A k . At tree-level these phases appear only 
in certain combinations, however, and exploiting tadpole conditions we are left with only one 
independent phase combination. Explicit CP violation in the Higgs sector leads to potentially large 
corrections to the electric dipole moments (EDMs). The non-observation of EDMs for thallium, 
neutron and mercury [17] severely constrains the CP-violating phases. However, as the phase 
combinations occurring in the EDMs can be different from the ones inducing Higgs mixing, the 
phases can be chosen such that the contributions to the EDMs are small while the phases important 
for the Higgs sector can still be sizeable |18| . This provides additional CP violation necessary for 
electroweak baryogenesis |19] . The explicit tree- level CP violation induces scalar-pseudoscalar 
mixings between the doublet fields H u ^ and the singlet field 5, but not between the scalar and 
pseudoscalar components of the Higgs doublets H u ^- The latter is realized in scenarios where 
explicit CP violation in the Higgs sector is induced through radiative corrections. Radiative CP 
violation furthermore allows for a moderate amount of CP violation which is still in agreement 
with the bounds from the EDMs |20|. In this respect, the CP phases which play a role are <pA t % PA,, 
from the trilinear couplings At,Ab. They are involved in the dominant corrections from the third 
generation squark loops. Phases from third generation Yukawa couplings on the other hand can 
be reabsorbed by redefinitions of the quark fields when neglecting generation mixing. Further 
sources for radiative CP violation stem from the gaugino sector where the soft SUSY breaking 
mass parameters M\,M2 and M3 are in general complex. One of the two parameters M\ and M2 
can be chosen real by applying an i?-symmetry transformation. The gluino mass parameter M3 
and hence its phase enters only at the two-loop level. 

Radiatively induced CP-violating effects from the third generation squark sector have been 
considered in the effective potential approach at one-loop level in Refs. [21]. One-loop contributions 



from the charged particle loops have been taken into account by Refs. 22 , also in the effective 



potential approach. In Ref. 23 the third generation (s) quark and gauge contributions are included 
in the one-loop effective potential. The full one-loop and logarithmically enhanced two-loop effects 
in the renormalization-group improved approach have been included in [24] . In order to properly 
interpret the results from the experiments and distinguish the various Higgs sectors from each other, 
a precise knowledge of the Higgs boson masses and couplings at the highest possible accuracy, 
including higher-order corrections, is indispensable. In this paper we consider the full one-loop 
corrections to the Higgs boson masses in the CP-violating NMSSM in the Feynman diagrammatic 
approach^] We allow for explicit CP violation at tree-level by including non-vanishing CP phases 
for the Higgs doublets and singlet, as well as for A, k, A\ and A K , which effectively reduce to 
one physical CP phase combination at tree-level when the tadpole conditions are exploited. We 



2 Higher-order corrections to the Higgs boson masses in the real NMSSM can be found in Refs. 



2 r > 
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furthermore allow for radiatively induced CP violation stemming from the stop)^] by choosing a 
non- vanishing CP phase for At. The Higgs sector as well as the neutralino and chargino sector 
will be used to determine the counterterms. The renormalization is performed in a mixed scheme 
which combines DR conditions for the parameters not directly related to physical observables with 
on-shell (OS) conditions for the physical input values. For the choice of our parameter sets the 



recent constraints from the Higgs boson searches at LEP 27 , Tevatron |28| and LHC 111 2] are 



taken into account. The inclusion of CP violation affects the Higgs phenomenology and hence the 
validity of possible scenarios. Our results will therefore help to clarify the question what kind of 
Higgs sector may be realized in nature for a SM-like Higgs boson with mass around 125 GeV, should 
the tantalizing hints of the LHC experiments be confirmed by the discovery of the Higgs boson at 
the 5<T level once a sufficient amount of data is accumulated. 

The organization of the paper is as follows. In Sect. [2] the parameters of the complex NMSSM 
will be introduced. Section [3] presents the details of our calculation. After introduction of the 
complex tree-level Higgs sector in Sect. 3.1 the set of input parameters is given in Sect. 3.2. The 
renormalization conditions and determination of the Higgs masses as well as the mixing angles are 
discussed in Sects. 3^5 and |3.4[ respectively. In Sect. [4] we present our numerical analysis, discuss the 
influence of different renormalization schemes as well as the consequences of one-loop corrections 
in the complex NMSSM for the Higgs boson phenomenology and results at the LHC. We terminate 
with the conclusions in Sect. [H 



2 Complex parameters in the NMSSM 

The Lagrangian of the complex NMSSM can be divided into an MSSM part which is adopted 
from the MSSM Lagrangian and an additional NMSSM part. The latter contains (apart from the 
phases of the Higgs doublets and singlets) four additional complex parameters. Two of them are 
the coupling k of the self-interaction of the new singlet superfield S and the coupling A for the 
interaction of S with the Higgs doublet superfields H u and H d (H u and H d couple to the up- and 
down-type quark superfields, respectively). They are introduced via the extension of the MSSM 
superpotential Wmssm> 

Wnmssm = Wmssm — eab^SH^Hu + £ K S 3 > (1) 

with ei2 = e 12 = 1. The MSSM part of the superpotential is given by 

V^mssm = -e ab {y u H a a Q b U c - y d H a d Q h D c - y e H a d L b E c ) , (2) 

where Q and L are the left-handed quark and lepton superfield doublets and U, D and E are the 
right-handed up-type, down-type and electron-type superfield singlets, respectively. The super- 
script c denotes charge conjugation. Colour and generation indices have been omitted. The quark 
and lepton Yukawa couplings are given by y d , y u and y e . They are in general complex. However, 
when neglecting generation mixing, as we assume in this paper, the phases of the Yukawa couplings 
can be reabsorbed by redefining the quark fields i.e. these phases can be chosen arbitrarily without 
changing the physical meaning [29] . 

3 For the low values of tan /3 applied in our numerical analysis, the CP-violating effects from the sbottom sector 
are marginal. 
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The soft SUSY breaking Lagrangian in the NMSSM is also extended with respect to the MSSM, 
^nmssm = £mssm " ™ 2 S \S\ 2 + {e ah A x \SH a d H b u - ±A k kS 3 + h.c.) , (3) 

containing two further complex parameters specific to the NMSSM, the soft SUSY breaking trilinear 
couplings A\ and A K . The MSSM part is given bjQ 

/•soft 2 \Tj \2 2 \tj \2 2 I /Si 2 2 \ ~ \2 2 | J |2 2 | f |2 2 |~ |2 

Aissm = - m H d \ H d\ -m Hu \H u \ - m Q \Q\ - m^un] - m D \d R \ - m L \L\ - m E \e R \ 
+ e ab {y u A u H a a Q b u R - y d A d H2Q b d* R - y e A e H a d Q b e* R + h.c.) 

- X - {M\BB + M 2 WiWi + M 3 GG + h.c) . (4) 

The soft SUSY breaking trilinear couplings A u , A d , A e of the up-type, down-type and charged 
lepton-type sfermiona^J respectively, which are already present in the MSSM, are in general complex. 
However, the soft SUSY breaking mass parameters of the scalar fields, m 2 x (X = S, H d , H u , Q, 
U, D, L, E), are real. The SM-type and SUSY fields forming a superfield (denoted with a hat) 
are represented by a letter without and with a tilde, respectively: Q, L and u R , d R , e R are the 
superpartner fields corresponding to the left- and right-handed quark and lepton fields. In general, 
also the soft SUSY breaking mass parameters of the gauginos, Mi, M2 and M3, are complex 
where the gaugino fields are denoted by B, Wi {i = 1,2,3) and G for the bino, the winos and the 
gluinos corresponding to the weak hypercharge U(l), the weak isospin SU(2) and the colour SU(3) 
symmetry. The /^-symmetry can then be exploited to choose either Mi or M2 to be real. The 
kinetic and gauge interaction part of the NMSSM Lagrangian finally do not contain any complex 
parameters. 

Expressing the Higgs boson fields as an expansion about the vacuum expectation values, two 
further phases appear, 

~-(^ Vd + hd + m A H -e^( K \ S- e —(v+h+ia) (5) 

Hd ~{ h- ) ' Hu ~ e {^(vu + h u + ia u )) ' S ~ ^2 {Vs + hs + ta '>- (5) 

The phases (p u and (p s describe the phase differences between the three vacuum expectation values 
{H®), (H®) and (S). In case of vanishing phases, <p u = cp s = 0, the fields hi and en with i = d,u,s 
correspond to the CP-even and CP-odd part of the neutral entries of H u , H d and S. The charged 
components are denoted by hf (i = d,u). 

Exploiting that the phases of the Yukawa couplings can be chosen arbitrarily, the phase of the 
up- type coupling is set to (p Vn = — (p u while the down- type and the charged lepton-type ones are 
assumed to be real. This choice ensures that the quark and lepton mass terms yield real masses 
without any further phase transformation of the corresponding fields. 

In the following renormalization procedure we will make use of the chargino and neutralino 
sectors, therefore they are introduced briefly here. The fermionic superpartners of the neutral Higgs 



4 Here the indices of the soft SUSY breaking masses, Q (L) stand for the left-handed doublet of the three quark 
(lepton) generations and U,D,E are the indices for the right-handed up- type and down- type fermions and charged 
leptons, respectively. In the trilinear coupling parameters the indices u, d, e represent the up-type and down-type 
fermions and charged leptons. 

5 We neglect generation mixings so that we have nine complex numbers A u , Ad, A e instead of three complex 
3x3 matrices. 
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bosons and colourless gauge bosons are H% and H® for the neutral components of the Higgs doublets, 
S for the Higgs singlet, the bino B and the neutral component W3 of the winos. After electroweak 
symmetry breaking these fields mix, and in the Weyl spinor basis ip° = [B,Ws, H%, H®, S) T the 
neutralino mass matrix Mjv can be written as 

M N = 

( M l -c p M z se w Mzapa^e-*^ 

M 2 CfiM w -Mwspe-*** 

-Cf3M z se w cpM w 



_\v a e vp a ^M w spse w Xe i ^ 
V2M w CfjSg w X 



(6) 



\ V2M w spse w \e^ y/2M w cpse w \ ^ ilps j 

\ ' e e V s / 



where Mw and Mz are the W and Z boson masses, respectively. The angle (3 is defined via 
the ratio of the vacuum expectation values of the two Higgs doublets, tan/3 = v u /vd, dw denotes 
the electroweak mixing angle and e is the electric charge. From here on the short hand notation 
c x = cosx, s x = sinx and t x = tanx is used. 

The neutralino mass matrix Mjy is complejj^] but symmetric and can be diagonalized with the 
help of the 5 x 5 matrix M, yielding diag (m^o , m^o , m^o , m^o , m^o ) = J\f* MjyN^ , where the absolute 

ai A 2 A3 A4 A.5 y 

mass values are ordered as Im^ol < ... < Im^ol. The neutralino mass eigenstates x?> expressed as a 
Majorana spinor, can then be obtained by 

X° = (X°,X°) T with Xi=Mirtp i,j = l,-..,5. (7) 

The fermionic superpartners of the charged Higgs and gauge bosons are given in terms of the 
Weyl spinors , , W\ and W2 where the latter two can be reexpressed as = (Wi^iW^) / V%- 
Arranging these Weyl spinors as 

^ = IX- 1 > V£ = | l + ) (8) 
leads to mass terms of the form, (tp^) T + h.c, with the chargino mass matrix 





The chargino mass matrix can be diagonalized with the help of two unitary 2x2 matrices, U and 
V, yielding diag(m~±, m~±) = U*McV* with m~± < m~±. The left-handed and the right-handed 

Al A2 Al A2 

part of the mass eigenstates are 

xi = v^+, r R = ur R , (10) 

respectively, with the mass eigenstates xf = (Xl-j Xr ) T > « = 1, 2, written as a Dirac spinor. 
6 Note, that in general the parameters A, k, Mi and M2 are complex. 
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3 The Higgs Boson Sector in the Complex NMSSM 



3.1 The Higgs Boson Sector at Tree-Level 



To ensure the minimum of the Higgs potential Vniggs at non-vanishing vacuum expectation values 
v ui v di v s the terms linear in the Higgs boson fields have to vanish according to 



|Min. 



for (j) = h d ,h u ,h s ,a d ,a u ,a s . 



(11) 



At tree-level, these tadpole parameters t$ are given bjj^] 

tK = K. - % - ^<^W + N f cv.) + |A| 



2 



1.4 



«2is|A|(— ^-Cp. + M^<vJ - |A| w s 



2MiW a 



1 



W VV i I |2 3 , 1 u II I 2 

' \ K \ V s ' ^r^\ /i K\\K\V s C { p z , 



1 



t 



-t 



1 



2My V s 2 dw S2p 1 , 

ta s = ^ 1^1(^=1^15^ + \K\v s S(p y ) — -^\A k \\k\v s Sip z , 



where we have introduced a short hand notation for the following phase combinations 

fx = fA x + fX + fs + fu , 
fy = fn ~ fX + 2<Ps ~ fu , 
fz = fA n + fn + 3fs ■ 



(12) 
(13) 

(14) 

(15) 
(16) 

(17) 



(18) 
(19) 
(20) 



In the expressions for the tadpole parameters some of the original parameters have already been 
replaced in favour of the parameters on which we will impose our renormalization conditions, as 



described in detail in Sects. 3.2 and 3.3 Thus the vacuum expectation values v u , v d and the U(l) 
and SU(2) gauge couplings g' and g have been replaced by tan/3 = v u /v d , the gauge boson masses 



M\y and Mz and the electric charge e (according to Eqs. (124) and (125) in Appendix |A[). This 
replacement has also been applied in the expressions of the mass matrices given below. 



It should be noted that the Eqs. (15) and (17) can have zero, one or two solutions for f x ,f z G 



[—ir, 7r) depending on the values of the parameters. If no solution is found there is no minimum 
of the Higgs potential at the corresponding set of values v d ,v u , v s and thus this parameter point is 
discarded. The single solutions yield one of the two values f x ,fz = ±vr/2. Assuming there exist 



two solutions of Eq. (15) then if with > solves this equation then also it — tp~. is a solution 



and similarly if tp~ with ip%. < is a solution then — (it — (p~,) is the second solution, analogously 
for <pf . 



7 The complex parameters are expressed in terms of their absolute value and a complex phase, i.e. for example 
A = |A|e i¥, \ 
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The terms of the Higgs potential which are bilinear in the neutral Higgs boson fields contribute 
to the corresponding 6x6 Higgs boson mass matrix in the basis of (j) = (hd, h u , h s ,a,d, a u , a s ) T 
which can be expressed in terms of three 3x3 matrices Mhh, M aa and M^ a 



Ma, 



M hh M ha 
Ml M aa/ 



(21) 



where Mhh and M aa are symmetric matrices. 

The entries of Mhh describing the mixing of the CP-even components of the Higgs doublet and 
singlet fields read 



M h d h d = Mfc| + -\\\v s tp(V2\A x \c Vx + \K\v s c Vy ) , 

11 M 2 s 2 

M hd h u = — 2 M\s 2fi - -\X\v s (V2\A x \c Vx + \K\v s c Vy ) + 2\X\ 2 ^^ 

M huhu = M 2 z s} + -\X\^(y/2\A x \ Ctpx + \k\v sCiPv ) , 



S2/3 



M hd h B 



2 A -cqv s - A 

e 

M huh3 =2\X\ 2 ^^s p v s -\X\ 



fx 

M W S9y 

e 

M w s e , 



-s l3 (V2\A x \c tpx + 2\k\v s c. 



-cp{y/2\A\\c Vil! + 2\k\v s c (Pv ) 



M hshs = 2\k\ 2 v 2 + ^=\k\\A k \c^ + V2\A X \ 



M w4 

j W tt\ 

e 2 v s 



(22) 

(23) 
(24) 

(25) 
(26) 

(27) 



Note, that here the tadpole conditions Eqs. (12)— (14) together with Eq. (11) have already been 
applied to eliminate , m2 H u an d Trig. Exploiting additionally Eqs. ( 15 ) and (17) we can eliminate 
c Vx and c Vz through 



±4/1 



k\ 2 v 2 „ 
— — -s 2 



2\A X \ 2 



±4/1- 18 



eWv 2 



(28) 



(29) 



The two signs correspond to the two possible solutions of Eqs. (15) and (17) (as explained before). 
Choosing either solution will define the sign of the cosine. In our numerical evaluation we will treat 
the sign as a further input. 

The mixing of the CP-odd components of the Higgs doublet and singlet fields is characterized 
by the matrix M aa which has the following entries, 



M, 



X\(V2\A x \c Vx 



M. 



a d a u 



M, 



1 
2 

c <Px + \ K \v s c Vy )v s 
-\\\(V2\A x \ Cipx + \K\v s c Vy )j-^ 
,M w s e 



M adas = \X\^J^ Sp {V2\A x \c^ - 2\k\v s c Vv ) , 



(30) 
(31) 
(32) 

(33) 
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M au a 3 = |A| —^cp 



{y/2\A x \t 



2\k\v.c, 



ipyl > 



M 2 s 2 

M asas = |A|(\/2|A A | C ^ +4\K\v sC(f v ' " H,r 



e 2 v. 



S2/3 - 31^1^1-^1^ , 



(34) 
(35) 



where again Eq. (11) together with the Eqs. (12)-(14) have been applied and Eqs. (11), (15) and 
(17) can be used to replace c lfx and c ipz . The matrix M^ a governs the mixing between the CP-even 
and the CP-odd components of the Higgs doublet and singlet fields, 

/ 3v s S/3 \ 

3v s c, P^MIAI^ 

A M w s I ' 
\-VsSj3 -V s C/3 -4s 2 /3 d - 



M, 



ha 



(36) 



In case of ip y = n y ir, n y 6 Z, the entries of Mh a vanish and hence, in that case, there is no CP 
violation at tree-level in the NMSSM Higgs sector; after transformation to the mass eigenstates we 
are left with three purely CP-even and two purely CP-odd Higgs bosons. 

The transformation into mass eigenstates can be performed in two steps. In our approach, first, 
the Goldstone boson field is extracted by applying the 6x6 rotation matrbfj7£ G , 



where $ = (hd,h u , h s , A,a s ,G) T . The resulting mass matrix, 

M$$ = K g M h K gT , 
can finally be diagonalized with the help of the matrix 1Z leading to 



diag((MW) 2 ,...,(Mg) 2 ,0 



--:V H 



(37) 
(38) 
(39) 



with the mass values being ordered as < ... < and the superscript (0) denoting the 

tree-level values of the masses. The corresponding mass eigenstates are obtained as 



H; 



The mass matrix M h + h - of the charged entries of the Higgs doublet fields, 



is explicitly given as 



M h+h - 



1 tt 



(h + d ,h+)M h+h -{h- d ,K) T 



MyyS2f3 + IAIv^V^AaIc^ + \k\v 8 C 



<fy, 



M 2 S 2 

2|A| 2 =^ 



S2/3 



(40) 
(41) 

(42) 



where again the Eqs. (11 )— ( 14 ) have already been applied. Diagonalizing this mass matrix with 
the help of a rotation matrix with the angle /3 C , where /3 C = (3 at tree-level, yields the mass of the 
physical charged Higgs boson, 

Mir si 



\X\Vs 

S2/3 



(V2\A X \, 



■fx 



+ \k\v s c, 



l w°e w 



(43) 



and a mass of zero for the charged Goldstone boson. 
8 The explicit form of 1Z G can be found in Appendix B 
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3.2 Set of Input Parameters for the Higgs Boson Sector 

To summarize, the original parameters entering the Higgs potential and thereby also the Higgs 
mass matrix are 

m H d i m H v ,i m s><PA K ,<PA x , \A\\,g,g',v u ,v d ,v s ,(p s ,ip u , \\\,<p\, \k\,<p k , \A K \ . (44) 

Instead of using this set of original parameters it is convenient to convert it to a set of parameters 
which offer an intuitive interpretation. This is especially true for the parameters which can be 
replaced by gauge boson masses squared as they are measurable quantities. We have chosen the 
set, 



t hd ,t hu ,t hs ,t ad ,ta s ,Mg ± ,M^,Mz,e,t&nP,v s ,(p s ,(p u , |A|,c/? A , N,<P*, \A K \ (45) 
s v ' v -v- ' 

on-shell DR 

where the first part of the parameters are directly related to "physical" quantities^] and will be 
defined via on-shell conditions while the remaining parameters are understood as DR parameters 
(see Sect. [3lT2[ ). The transformation rules for going from set Eq. (44) to set Eq. (45) are given in 
Appendix |Aj 

3.3 The Higgs Boson Sector at One-Loop Level 

At one-loop level, the Higgs boson sector and the corresponding relations between parameters of the 
theory and physical quantities are changed by radiative corrections. In particular, the Higgs boson 



mass matrix receives contributions from the renormalized self-energies 10 Sy(p 2 ) at an external 
momentum squared p 2 , 

± i3 (p 2 ) = Z i:j (p 2 ) + ^p 2 \bZ^ + 5z] - \ \dZ^V H + v\jSz] . - [KSMttR^y , (46) 



2 



with i,j = 1, ... ,6 and Hq = G the Goldstone boson. The unrenormalized self-energies Sy are 
obtained by taking into account all possible contributions to the Higgs boson self-energy, including 
the ones from fermion, gauge boson, Goldstone boson, Higgs boson, chargino, neutralino, sfermion 
and ghost loops. 

The wave function renormalization matrix 5Z in the basis of the Higgs boson mass eigenstates 
is derived via rotation from the corresponding matrix 5Z§ in the basis of the Higgs boson states 

*, 

bZ = U5Z^n ] . (47) 

The Higgs boson fields $ are renormalized by replacing the fields by renormalized ones and a 
renormalization factor. This can be expressed as, valid up to one-loop order, with the field renor- 
malization constant 5Z§ as 

$->(l+^Z*)$, (48) 

9 Whether the tadpole parameters can be called physical quantities is debatable but certainly their introduc- 
tion is motivated by physical interpretation. Therefore, in a slight abuse of the language, we are also calling the 
renormalization conditions for the tadpole parameters on-shell. 

10 In general, we call E and E renormalized and unrenormalized self-energy, respectively. 
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where 5Z§ = Tip bZ^TiP^ . The field renormalization constant 5Z^ of the interaction eigenstates 
is a diagonal matrix 



5Z<j> = d\a,g(5Z Hu) 5Z Hd , 5Z s ,SZ Hd , SZ Hu , 5Z S ) . 



(49) 



The explicit definitions and expressions for SZh u , SZjj, and 5Z$ are given in Sect. 3.3.1 



The matrix 5A4$$ denotes the counterterm matrix in the basis of the Higgs boson states <J> 
which has been introduced within the renormalization procedure by replacing the parameters in 
the mass matrix as given in Appendix |A| in Eqs. (118)-(125) by their renormalized ones and 
corresponding counterterms and expanding about the renormalized parameters. The part linear in 
the counterterms forms the mass matrix counterterm. The specific definitions of the parameters 
and the determination of the counterterms are discussed in Sect. I3.3.2L 



3.3.1 Higgs boson field renormalization 



The field renormalization constants introduced in Eq. (49) are defined in the DR scheme. The 
precise expressions for 5Zjj d , SZfj u and 5Z$ are determined via the following system of equations 



SZ Hi {\Tlix\ 2 + \n H sm/3 + K i6 cos/3| 2 ) + 6Z Hu {\K i2 \ 2 + |ft i4 cos/3 - 7^ 6 sin/3|' 



+ 8Z s {\K i3 1 2 + \TZ 



where 



i~> 



V' 



div 



with 



1,2,3 , 



dp 2 



p 2 = ( M (0))2 



(50) 



(51) 



The subscript 'div' denotes that only the divergent parts proportional to A are taken into account 
with A = 2/(4 — Z?) — 7£ + In 4-7T and being the Euler constant and D the number of dimensions. 



The pole of A for D = 4 characterizes the divergences. Solving this system of equations Eq. (50) 
yields 



with 



sz Hd 

5Z S 



m 

Til 

r 



[{r 23 r 32 - r 22 r 33 )T,' u + (ri 2 r 33 - ri 3 r 32 )£' 22 + (r 13 r 22 - ri 2 r 23 )£j 



33Jdiv ' 



[(r 2 ir 33 - r 23 r 31 )T,' n + (ri 3 r 31 - rnr 33 )S 22 + (r u r 23 - r 13 r 2 i)T!., 



33Jdiv ' 



[{r 22 r 31 - r 2 ir 32 )T,' u + (rnr 32 - r 12 r 3 i)E' 22 + (r 12 r 21 - rnr 22 )£j 



33Jdiv ' 



(|^i| 2 + |^i4 smp + TZiQ cos/3| 2 ) 
(\U i2 \ 2 + \K i4 , cos/3 - Tlx sin/3| 2 ) 

(\n i3 \ 2 + \n i5 \ 2 ) , 

mr 22 r 33 + n 2 r 23 r 3 i + ri 3 r 32 r 2 i 



^11^23^32 - ^13^22^31 ~ ^12^21^ . 



(52) 
(53) 
(54) 



(55) 
(56) 
(57) 
(58) 



It should be noted that TZte = for i ^ 6 and hence, in Eqs. (55) and (56) terms proportional to 



TZiQ vanish for the values i = 1, 2, 3 needed in Eqs. (52)-(54). 
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3.3.2 Parameter renormalization 

The parameter renormalization is performed by replacing the parameters by the renormalized ones 
and the corresponding counterterms, 

t<j> -> tfj, + 5t^ with (p = {h d ,h u ,h s ,a d ,a s } , (59) 

M 2 H± -> M^± + 5M 2 H± , M^^M^ + 5M^ , M 2 Z -> M| + <5Mf , (60) 

e -»• (1 + 5Z e )e , (61) 

tan /3 — )• tan (3 + 5 tan /3 , i> s — > v s + 5v s , (62) 

->■ fs + Sifi s , -)• <£ u + <5</? u , (63) 

A -> A + 5A = A + e iv=,x (5| A| + i\5tp x , k ^ k + 5k = k + e ilfK 5\K\ + iK5tp K , (64) 

\A K \^\A K \ + 5\A K \ . (65) 

In the case of complex parameters the complex counterterms can be understood in terms of two 
real counterterms, one for the absolute value and one for the phase, as in Eq. (64) for 5X and 5k. 

As we make use of the chargino and the neutralino sector for the determination of the coun- 
terterms 5v s , 5(p s , 5X, 5k and 5ip u we also need to renormalize the gaugino mass parameters M\ 
and M 2 , 

Mi ->• Mi + SMi = Mi + e^ M i 6\Mi | + iM x 5^ Mx , (66) 
M 2 -»• M 2 + £M 2 = M 2 + e^ A/ 2(5|M 2 | + iM 2 5tp M2 ■ (67) 

To keep the relations as general as possible we do not make use of the i?-symmetry relations here 
and keep both gaugino mass parameters complex. 

In the following, we list all the renormalization conditions and counterterms. The renorma- 
lization scheme applied here is a generalization of the "mixed scheme" of Ref. 26 for complex 



parameters - we will be brief on the conditions that can be directly taken from Ref. 26 

(i-v) Tadpole parameters: 

The renormalization conditions for the tadpole parameters are chosen such that the linear 
terms of the Higgs boson fields in the Higgs potential also vanish at one-loop level, 

5t <t> = T 4> with 4> = h d ,h u ,h s ,a d ,a s , (68) 

where denotes the contribution of the irreducible one-loop tadpole diagrams. 

- viii) Masses of the gauge bosons and the charged Higgs boson: 

The masses of the gauge bosons and of the charged Higgs boson are determined via on-shell 
conditions requiring that the mass parameters squared correspond to the pole masses squared 
leading to 

5M 2 V = ReY^ w (M^) , 5M 2 = ReZ T zz (M 2 ) , 5M 2 H± = ReX HTH± (M 2 H± ) , (69) 

where and ^ zz are the transverse parts of the unrenormalized W boson and Z boson 

self-energy, respectively, while Ti H ^ H ± denotes the unrenormalized self-energy of the charged 
Higgs boson. Re takes only the real part of the scalar loop functions but keeps the complex 
structure of the parameters. 
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(ix) Electric charge: 

The electric charge is fixed via the ee'j vertex in such a way that this vertex does not receive 
any corrections at the one-loop level in the Thomson limit, i.e. for zero momentum transfer. 
This yields (cf. Ref. [30] up to a different sign convention in the second term) 

5Z e = ±s£(0) + -^££(0) , (70) 

z c 8w lvl Z 

with E^ and T^ z being the transverse part of the unrenormalized photon self-energy and 
the unrenormalized mixing of photon and Z boson, respectively. 

(x) Ratio of the vacuum expectation values tan/3: 

The ratio of the vacuum expectation values tan j3 is defined as a DR parameter and the 
counterterm is given by |31l|32| 



5 tan/3 = [5Z Hu - 5Z Hd ] | div . (71) 

(xi,xii) Vacuum expectation value v s and phase ip s : 

The singlet vacuum expectation value v s and the phase ip 8 are determined as DR parameters. 
For the derivation of the corresponding counterterms, we start out from the on-shell conditions 
for the chargino masses, 

ReS x +(p) xt(p)\ p 2= m \ =0, i = 1, 2, (72) 

where E + are the renormalized chargino self-energies. Applying the decomposition of the 
fermionic self-energy 

% 3 {p 2 ) = ^( P 2 )V L +i>±*{p 2 )V R + tff(p 2 )P L + t^(p 2 )P R (73) 

with Vl,r = (l : F75)/2 being the left- and right-handed projectors, leads to the finite relations 

[m-± (ReE^ ± (p 2 ) + ReE« ± (p 2 )) + ReE^ (p 2 ) + ReE^ (p 2 )]^ = , (74) 

[m-± (ReS^ ± (p 2 ) - ReS J ± (p 2 )) - ReE^ (p 2 ) + ReSjl (p 2 )] . . = , (75) 

which can be exploited for the determination of the counterterms. Using the expressions 
given in Eqs. ( |143| )-( [T46 ) in Appendix [C] for the renormalized self-energies yields 

Re([/W c yt)| d . v = 1 [ m . f (^ + (p 2 ) + Sj + (p 2 )) + £^(p 2 ) + S^(p 2 )] ..| div 

=: Re 5m + , (76) 



lra(U*5M c V^)\ div = [S^(p 2 ) - E^(p 2 ) + imMU* \^5Z C R U T + VIm5Z<£v% 



Idiv 



: Im 5m + . (77) 

A r ' 



The imaginary parts of the field renormalization constants have been set to zero, hence 
ImJZg = Tm5Z$ = 0. With 

/ 5M 2 y/^-^[S(M W 8p)-is p M w S<p u ]\ 
6M ° = \V2 5{c p M w ) ^ [X8v s + iXvJtps + v s 5\] ) ' 
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solving Eqs. (|76|) and <\77\) for Sv s + iv s Sip s and SMi 11 we obtain 



Uh = \U u \ 2 -\V 12 \^ VllUu6m ^ " V2lU2l6m xt 2 

- U n U^SM C21 \ div - VnV&SMcnUv] , 

\pl\*e~ iLPa 

5v s + iv s 8ips = m2 ,,, r l2 2 [ - V12C/12 5m + + V22^22 5m + 

|A| 2 (|C/n| 2 - |Vi 2 | 2 ) L Xl1 * 22 

+ ^C/ia 5M Cl2 | div + VftVia 5M C21 | div ] - u s A* 



5A 



(79) 



(80) 



It should be noted that 6<p u contained in 5Mc 12 as well as 5X have not been defined yet. We 
need additional conditions given by Eqs. (98)-(100). Together with Eqs. (79) and (80) they 
form a system of linear equations that can be easily solved but leads to lengthy expressions. 

(xiii - xvii) Couplings A and k and the phase ip u : 

The phase ip u as well as A and k are also defined as DR parameters. On-shell conditions for 
the neutralino masses 



12 



=m n 



0, i = 1 

are exploited to derive the following finite relations^' 



j • • • 5 ^ j 



(81) 



[m^o (ReSj (p 2 ) + ReEj (/)) + Ret^(p 2 ) + ReE^V)]^ = , 
[m^o(Re^ (p 2 ) - Re£* (p 2 )) - Re± L J( P 2 ) + ReSjf(p 2 )]^ = , 



(82) 
(83) 



with i = 1, . .. ,4. Applying the Eqs. ( 150 H 153) in Appendix |C| leads to 

Re(M*6M N M^ii\ div = [m^ (p 2 ) + \{^(p 2 ) + (p 2 ))] ..\ div =: Re5m x o , (84) 



lm(Af*8M N Af%\ Aiv = - [Sjf(p 2 ) - E£oV) + 2im^(A/'*Im5Z JV A/'t)]. i 
=: Im5m v o , 



Idiv 



(85) 



i = 1, ... ,4 where already £ |jj = Tj^ \h has been used which is true due to the Majorana 
character of the neutralinos. The imaginary part of 5Z N has been set to zero, Im^Z^ = 0. 
The elements of the mass matrix counterterm 5Mn are derived as 



5M Nll = 5Mi 
SM N22 = 5M 2 



(86) 
(87) 



11 Evon though M2 does not enter the Higgs boson sector at tree-level, it has to be dealt with due to its entanglement 
in Eq. (781. 

12 In terms of an on-shell scheme, if both chargino masses were defined on-shell only three of the neutralino masses 
could be chosen independently. Although Eq. (81 I leads to two independent equations for each i, for i = 4 it is only 
partly used to fix one still undefined phase. 

13 It should be noted that, in general, both Eqs. (821 and (83 1 for i = 4 only hold for the divergent part simultane- 
ously as all the parameters are already fixed by other conditions. 
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5M N55 
8M Nl3 

5M Nl4 
6M N23 
5M N24 

5M N:j4 



5M N35 



8M Ni5 
SM N33 



V^e^" [v s 5k + k5v s + iK,v s 5(p s ] , 

-M z se w cl S p 6 tan (3 + °\ [5M } 

2s dw M z 

e-^«[5(s /3 M z s ew ) - ispM z s 0w 5^ u 
S(cpM w ) , 

-e-**»[8{spM w ) ~ ispM w 5ip u ] , 



w ~ SM 2 Z ] , 



1 



e tLPa [v s SX + X5v s + i\v s 5ip s ] 



[x 5{s p M w s dw ) _ s p M w s dw ^ _ ^ _ 



2 A 



8{cpM w s 6w ) cpMwse 



(5X-XSZ e 



+ 

e e 
5Miv 44 = 5M Wl2 = 5M Nl5 = SM N25 = . 



(88) 
(89) 

(90) 
(91) 
(92) 

(93) 
(94) 

(95) 
(96) 



As the neutralino mass matrix is symmetric, this also holds for the counterterm mass matrix 
and therefore SM^.. = 5M^. V Rewriting Eqs. (84) and (85) explicitly and solving for 5 Mi 
results in the following set of equations, 



5M 



AAA 



IMil 



«5m x?i - 2Af^ 3 5M Nl3 +Ml i 5M Nl4 \ - 2A/? 2 [A/? 3 <M% M + N^8M Nai 

- 2MU^uSM Ns4 +N? 5 5M N35 ] - 2M^Ml 5 5M N45 

- {M* l2 ) 2 5M 2 - (M^ 5 ) 2 5M N55 , 

2[a 2li 5M Nl4 + a 22i 8M N24 + 2a 23i 5M N34 + a 2 35<^v 35 + a 2i5 5M N45 ] 
+ a 222 5M 2 + a 255 5M N55 

= (M 21 ) 2 5m x0ii - (M^) 2 5m x o 2 - 2a 213 5M Nl3 - 2a 223 5M N23 , 
2[a3i 4 5MAr 14 + a 32i 5M N24 + 2a 334 5M Nz4 + a 335 5M N35 + a 345 (5MAr 45 ] 
+ a 322 SM 2 + a 355 5M Ns5 

= (A/" 3 *i) 2 5m x?i - {Mtx) 2 8m x ^ - 2a 313 5M Nl3 - 2a 323 5M N23 , 
Im|2[a 4 i 4 (5MAr 14 + a i2i 5M N24 + 2a i3i 5M N34 + a i35 5M N35 + a^ 5 5M N45 ] 

+ a^ 22 bM 2 + a455<->Miv S6 } 

2a423<5MAr 23 | , 



(97) 



(98) 



(99) 



= Im^X^fSm^ - (A/Ti) 2 £m x o 4 - 2a il3 5M N - 

where we have introduced the shorthand notation 

r* \2 * r* \r* t \ r* ^2 » r* » r* 



(100) 



a ijk = WYJWik ~ WiY^tk ■ 



(101) 



As stated above, taking Eqs. (98)-(100) together with Eqs. (79) and (80) leads to a system 
of equations with 4 complex and one real equation linear in the counterterms that has to be 
solved for 5M 2 , Sv s , 5ip s , 5X, 5k and 5ip u - 

(xviii) Absolute value of the singlet trilinear coupling \A K \: 

The absolute value of the singlet trilinear coupling \A K \ is determined as a DR parameter. 
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The corresponding counterterm is calculated using 

n i5 n j5 t ij (M asas ) = o, (102) 

which is equivalent to 

5M asas = ni S K j BEi j (M a ,a.) , (103) 
with 5M aaas depending on <5|j4 k |. Dropping the finite parts and solving for 6\A K \ yields 

s\A K \ = -— — z — — j—in^n^^Ma^-dh 



with 



\k\v s [3c Vz — 3t (Pz ^^\ *• 

3 3 i 

+ -^I^kI [ V sC v J\k\ + \K\CpJVs] + —\K\v s t lf Jf 2 \ (104) 
V 2 V 2 - 1 div 



r 2 ,,2 2 1 M ^ S L S i/3 M ^S^S 2 /3^c2 b ^ , _ 2 ^ 1 , t hs 



/l - [M H± - M w c Ap \ — j- 2 [t hu + tf}tp B t hd \ + 

+ \X\M^ S -^^[2\X\M^ v \ S2 p + 3\K\v 2 s c iPy ] , (105) 

gZ^jZ L gZ f r » J 

r2 „2 „_ . 

(106) 

The counterterms <J/i and 5/2 , which are functions of counterterms of the parameters defined 



3\\\M^si w s 2 f, 1 



as input in Eq. (45), are determined by replacing the parameters by renormalized ones plus 
corresponding counterterms and expanding about the parameters. It has to be taken into 
account that in the expressions A/3 = (3 — /3b only ft is treated within the renormalization 
procedure. The angle /3b = Pc = Pn 1S the mixing angle of the charged Higgs bosons and the 
angle extracting the Goldstone boson defined in Appendix |A} 

Following the approach above, it has been found that the counterterms S(p s , S(px, 5ip K , 5ip u , 



ScpMi and SipM 2 vanish. In that respect, it is interesting to note that Eqs. (77) and (85) allow for 
a certain freedom of choice; some potentially divergent parts can be moved into ImSZ 1 ^, ImdZ 1 ^ 
and Im^Z^ which do not appear in the calculation of the Higgs boson masses as the charginos and 
neutralinos only enter through internal lines in the Feynman diagrams. 

The derivation of the counterterms for v s , <p s , A, k and (p u presented above is not unique. In 
a second approach on-shell conditions for all the neutralino masses plus an additional condition 



from the chargino sector, Eq. (77) for i = 1, have been exploited to calculate the DR counterterms 
leading to the same result and providing a good cross-check. A further possibility is to determine 
the counterterms within the Higgs boson sector only. We have also done that but this does not 
test the calculation at the same level as using conditions from the chargino and neutralino sector. 

3.4 Loop Corrected Higgs Boson Masses and Mixing Matrix Elements 

The one-loop corrected scalar Higgs boson masses squared are extracted numerically as the zeroes 
of the determinant of the two-point vertex functions T, 

t{p 2 ) = i[l-p 2 -M lx ) with [M\ = (Mj®) 2 6 ij -t ij (j?) i,j = l,...,5, (107) 
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where Xjj(p 2 ) is given in Eq. (46). The superscript 11 denotes the one- loop order. 



Starting from Eq. (107) the Higgs masses at one-loop level can be obtained via an iterative 



procedure To calculate the one-loop mass of the n m Higgs boson the external momentum squared 
p 2 in the renormalized self-energies is set equal to the tree- level mass squared (p 2 = (M^) 2 ) 
in the first iteration step. Then, the mass matrix part of T, i.e. A4 , is diagonalized. The thus 
obtained n th eigenvalue is the first approximation of the squared one-loop mass. In the next 
iteration step p 2 is set equal to this value and once again the eigenvalues of A4 11 are calculated to 
yield the next approximation of the one-loop mass. This iteration procedure is repeated until a 
precision of 10 -9 is reached. All five Higgs boson masses are calculated this way. 



Note that in Eq. (107) the mixing with the Goldstone bosons is not taken into account but we 



have checked numerically that the effect is negligible. Furthermore, it was shown in Ref. 33 , that 
in the MSSM it is sufficient to include the mixing with the Goldstone boson, whereas the mixing 
with the longitudinal component of the Z boson does not have to be added explicitly. Taking into 
account the mixing of the Goldstone boson as well as the mixing of the Z boson leads to the same 
result as only including the Goldstone boson mixing. 

Due to the radiative corrections, not only the masses of the particles receive contributions but 
at the same time the tree-level mass eigenstates mix to form new one-loop mass eigenstates. In 
order to take this into account the Higgs mixing matrix 1Z which performs the rotation from the 
interaction eigenstates to the mass eigenstates has to be adjusted so that 

Hf = Kl& . (108) 

In the numerical analysis, for the simplicity of the notation we drop the superscript 11 again. If 
not explicitly mentioned one-loop corrections are included. 

The rotation of the tree-level to the one-loop mass eigenstates could be obtained by calculating 
finite wave function correction factors. The procedures to calculate these for a2x2or3x3 mass 
matrix are described in Ref. [13| and need to be extended for the 5x5 case. Another option 
is to apply the p 2 = approximation. After setting the momenta in A4 n to zero, the rotation 
matrix that relates the tree-level to the one-loop mass eigenstates can be defined as the matrix 
that diagonalizes A4 11 . The latter procedure has the advantage that the mixing matrix is unitary. 
The drawback is that it does not retain the full momentum dependence. For our numerical analysis 
we used the p 2 = approximation for the determination of the mixing matrix. But we checked 
numerically that the differences between both methods are negligible. 

4 Numerical Analysis 

The calculation of the one-loop corrected Higgs boson masses has been performed in two different 
calculations. While in one calculation the Feynman rules have been derived from the NMSSM 
Lagrangian and implemented in a FeynArts model file |34|, they have been obtained with the 



Mathematica package SARAH 35 in the second calculation and cross-checked against the first cal- 



culation. The self-energies and tadpoles have been evaluated with the help of FormCalc |36| in the 



14 This procedure is not strictly of one-loop order. It was shown, however, in Ref. |32| for the MSSM that this 
procedure gives much exacter values for the Higgs mass including implicitly higher order corrections than a strict 
treatment at one-loop level. 
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't Hooft-Feynman gauge. The divergent integrals, regularized in the constrained differential renor- 



malization scheme [37], have been computed numerically with LoopTools 36 . For the evaluation 
of the counterterms, numerical diagonalization of the one-loop corrected Higgs boson mass matrix 
and the determination of the mass eigenvalues finally two independent Mathematica programs have 
been written. 



We follow the SUSY Les Houches Accord (SLHA) 38 and compute the parameters Mjy and e 



of our input set defined in Eq. (45) from the SLHA pre-defined input values for the Fermi constant 
G F = 1.16637 • 1CT 5 GeV~ 2 , the Z boson mass M z = 91.187 GeV and the electroweak coupling 
a = 1/137. If not stated otherwise, we use the running DR top quark mass nit at a common scale 
Q = ^/m,Q 3 mt R - It is obtained from the top quark pole mass Mt = 173.2 GeV by taking the routines 
of NMSSMTools [39] . In the same way we obtain the running DR bottom quark mass starting from 
the SLHA input value mf,(m),) MS = 4.19 GeV. For the light quarks we chose m u = 2.5 MeV, 
m c = 1.27 GeV, m d = 4.95 MeV, m s = 101 MeV [40] and for the r mass m T = 1.777 GeV. 

In the following we exemplify the effects of complex phases in the one-loop corrections to the 
NMSSM Higgs boson masses in different scenarios. We require the scenarios to be compatible with 
the recent results of the LHC Higgs boson searches [lip] in the limit of the real NMSSM, i.e. for 



vanishing CP-violating phases 15 Our starting points are the NMSSM benchmark points presented 
in Ref. [8] which have been slightly modified for our analysis. With not too heavy stop masses and 
not too substantial mixing they avoid unnaturally large finetuning, and A and k have been chosen 
such that unitarity is not violated below the GUT scale. Furthermore, we paid attention to keep the 
effective fx < 200 GeV, with fx = Xv s /y/2, in order not to violate tree-level naturalness. For each 
scenario we verified that the non SM-like Higgs bosons are not excluded by the searches at LEP 27 
Tevatron 28 and LHC. This has been cross-checked by running the program HiggsBounds 4lJ 



it 



which needs the complex NMSSM Higgs couplings and branching ratios. The latter have been 



obtained by adapting the Fortran code HDECAY 46,47 to the complex NMSSM, in which we use 



the one-loop corrected Higgs boson masses and mixing matrix elements of our calculation. For 
the SM-like Higgs boson 

^SM-like of the real NMSSM we demand its mass to lie in the interval 
120 — 130 GeV. Furthermore, the total significance for Hf M ~ llke should not deviate by more than 
20% from the corresponding SM value. We roughly estimate the significance S to be given by 
S = N s /y/Nb, where N s denotes the number of signal events and 7V& the number of background 
events. Hence our criteria for a scenario to be compatible with present LHC searches are 

120 GeV < M^sM-iike < 130 GeV (109) 

5 NMSSM ( #SM-l* e) = S SM ( ^SM )±20% for M = M ^ . (HO) 

i 

In this case the scenario is estimated to be compatible with the present LHC searches taking into 
account experimental and theoretical uncertainties. We roughly approximate the total significance 
by adding in quadrature the significances of the various LHC Higgs search channels. We assume 
the number of background events to be the same both in the SM and the NMSSM case. For 
the calculation of the signal events we need the cross section values in the different channels, 



15 Note, that this is only an arbitrary choice which was made for practical reasons. We could as well have demanded 
a complex NMSSM scenario to be compatible with the recent LHC searches. 

16 The program NMSSMTools [39] also performs these checks. It can be used, however, only for the case of the real 
NMSSM. 
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which we obtain as follows. We first calculate the inclusive production cross section by adding the 
gluon fusion, weak boson fusion, Higgs-strahlung and ti Higgs production cross sections. Associated 
production with bb does not play a role here, as the tan j3 values we chose are rather low. The gluon 



fusion value at NNLO QCD is obtained with HIGLU 42 , which we have modified to the NMSSM 
case. Weak boson fusion and Higgs-strahlung at NLO QCD are computed with the programs VV2H 
and V2HV [43] by applying the modification factor due to the modified NMSSM Higgs coupling to 
gauge bosons compared to the SM case. Finally, the cross section value for ti Higgs production 
at NLO QCD [44] is obtained from the cross section values given at the LHC Higgs cross section 



working group webpage 45 by applying the appropriate factor taking into account the change of 
the NMSSM Higgs Yukawa coupling with respect to the SM coupling. Note that the NLO QCD 
corrections are not affected by changes due to the NMSSM Higgs sector and can therefore readily 
be taken over from the SM case. The cross sections in the WW, ZZ and 77 LHC search channels 
are obtained in the narrow width approximation by multiplication of the total cross section with 
the corresponding Higgs branching ratios into these final states. The branching ratios have been 



obtained from our modified Fortran code HDECAY 46 47 , adapted to the complex NMSSM. The 
thus obtained cross sections for the various channels can be used to calculate the number of signal 
eventsj^] The experiments take into account QCD corrections beyond NLO and also electroweak 
corrections. As these are not available for the NMSSM we cannot take them into account here. 
They are of the order of a few percent depending on the process. Furthermore, ATLAS and CMS 
exploit more final states and combine them in a sophisticated statistical procedure, while we have 
taken into account only the most prominent ones. Our approximation should therefore be viewed 
only as a rough estimate, good enough though to eliminate scenarios clearly excluded by the present 
LHC search results. 

In all investigated scenarios we have taken the input values at the scale Q = ^/mQ 3 mt R . In order 
to comply with the present LHC searches [48||49 , we have throughout taken the soft SUSY breaking 



mass parameters of the squarks of the first two generations equal to 1 TeV, and for simplicity also 
those of the sleptons. The corresponding trilinear couplings are taken to be 1 TeV. Furthermore, 
the right-handed soft SUSY breaking mass parameter of the sbottom sector is set equal to 1 TeV 
and its trilinear coupling close to 1 TeV, so that we have 

mu = mo = rriQ 12 = m E = m L = 1 TeV 
A x = 1 TeV (x = u, c, d, s, e, fi, r) 

A b « ITeV. (Ill) 

This leads to masses of ~ 1 TeV for the squarks of the first and second family, the sleptons and 
the heavier sbottom. Furthermore, all scenarios lead to the correct relic density in the limit of the 
real NMSSM, which has been checked with NMSSMTools which contains a link to MicrOMEGAs 501. 



17 Note, that the luminosity factor in the calculation of the number of events and also the number of background 
events drop out in the comparison of the NMSSM case to the SM case, so that we only need to calculate the quadratic 
sum of the cross sections in the different final states for the NMSSM and for the SM and compare them. 
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4.1 Scenario with a SM-like H$ 

The parameter set for this scenario is given by 

I A| = 0.72 , |k|=0.20, tan = 3, M^±= 629 GeV, |A K |=27GeV, = 198 GeV 

\A b \ = 963 GeV , \A t \ = 875 GeV , Mi = 145 GeV , M 2 = 200 GeV , M 3 = 600 GeV . 

(112) 

The slightly high values of A and k may require extra matter above the TeV scale |8| 18 We set 
all CP-violating phases to zero and subsequently turn on specific phases to study their respective 



influence. In this case, the signs of the tree- level CP-violating phases Eqs. (28), (29) are then 
chosen as 



sign cos tp x = +1 



sign cos tp z 



-1 . 



(113) 



Furthermore, the left- and right-handed soft SUSY breaking mass parameters in the stop sector 
are given by mg 3 = 490 GeV and mt R = 477 GeV. This leads to relatively light stop masses 
363 GeV and m; 



m 



fi 



616 GeV, still allowed by the experiments 



51 



52 
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In the calculation 



of the one-loop correction to the Higgs boson masses we have set the renormalization scale equal 
to 500 GeV, i.e. fj, ren = 500 GeV, if not stated otherwise. This scenario leads in the CP-conserving 
NMSSM to the one-loop corrected H3 being SM-like with a mass Mh 3 = 125 GeV compatible 
with present LHC searches. In the following we discuss for various complex phase choices the 
phenomenology of the three lightest Higgs bosons. The two heavier ones receive mass corrections 
of maximally 2 GeV leading to masses of ~ 642 GeV so that they are not excluded by present 
collider searches, with H4 being mostly CP-odd and H§ mostly CP-even. We therefore do not 
display their masses explicitly. 



4.1.1 CP violation at tree-level 



As we have seen in Sect. 3.1 a non-vanishing phase ip y , cf. Eq. (19), introduces CP violation at 
tree-level. Therefore CP-even and CP-odd Higgs mass eigenstates cannot be distinguished any 
more. A measure for CP violation concerning the state Hi (i = 1, 5) is instead provided by the 
quantity 



r'cv = (^i) 2 + {Ki2? + (Tt. 



(114) 



where IZij are the matrix elements of the mixing matrix which diagonalizes the Higgs boson mass 
matrix, cf. Eq. (39). A purely CP-even (CP-odd) mass eigenstate Hi corresponds to r l CP = 1 (0). 



We first investigate the effect of a non-vanishing phase 
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(115) 



The phases ^,,1^4,. are fixed by the tadpole conditions Eqs. (15), (17). 



1 Being above the TeV scale it is not expected to influence LHC phenomenology, apart from the indirect effect of 
allowing A to be a somewhat larger than allowed by the usual perturbativity requirement in the NMSSM with no 
extra matter. If instead one accepts more finetuning in the theory and allows for higher stop masses, a Higgs mass 
of the order of 125 GeV can be achieved for lower A values, cf. the discussion in Sect. 4.2 

19 Note, that light stop masses and small mixing reduce the amount of finetuning [8 . 

20 The choice of non- vanishing ip K allows to investigate mixing effects of the Higgs bosons while suppressing the 



phase relevant for the neutral electric dipole moment 23 
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In Fig. [T] (left) we show the tree-level and one-loop masses of the two lightest Higgs mass 
eigenstates H\ t 2 as a function of (p K , where (p K = corresponds to the real NMSSM. The phase is 
varied up to tt/8. Above this value it turns out that the phases ifA x and (p£ K cannot be chosen 
in such a way that the tadpole conditions are fulfilled. The tree-level and one-loop corrected mass 
of the SM-like H% is shown in Fig. [I] (right). Figure [2] displays, as a function of (p K , the amount 
of CP violation r l cp (left) and the amount of the CP-even singlet component]^] (7^3) 2 (right) for 
#1,2,3- Finally Fig. [3] shows their coupling squared to the V bosons (V = Z, W) normalized to the 
SM as a function of tp K . As expected, the masses exhibit already at tree- level a sensitivity to the 
CP-violating phase (p K . In particular for the SM-like H3 this dependence is more pronounced at 
one- loop level, changing its mass value by up to 9 GeV for ip K £ [0,7r/8]. The one- loop correction 
increases the mass by ~ 4 to 11 GeV depending on (p K with larger mass values for larger CP- 
violating phases. 

In the plots the grey areas are the parameter regions which are excluded due to the experimental 



constraints from LEP, Tevatron and LHC, and which have been obtained with HiggsBounds 22 This 



is the case for 0.074 tt < <p K < 0.099 tt and (p K > 0.112 tt. The dashed region excludes the parameter 



regions where the criteria stated in Eq. (110) of compatibility with the recent Higgs excess around 
125 GeV cannot be fulfilled any more, here for <p K > 0.0217T. The reason is that with increasing 
ip K the eigenstate H3 becomes more CP-odd and hence couples less to VV (V = Z, W) as can 
be inferred from Fig. [2] (left) and Fig. [3] so that the total cross section becomes smaller and the 
significance deviates by more than 20% from the SM significance of a SM Higgs boson with same 
mass. 

The one-loop corrections for the two lighter Higgs bosons H12 increase their masses by 6- 
15 GeV depending on (p K . The mass value M Hl r H2 \ decreases (increases) with rising (p K . In the 



21 The CP-odd singlet component is given by (IZis) 2 . 

22 The exclusion is due to the LEP constraint on Hi from the Higgs boson search in the Zbb final state stemming 
from a Higgs boson produced in Higgs-strahlung with subsequent decay into a fe-quark pair. 
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Figure 1: Left: Tree-level (dashed) and one-loop corrected (full) Higgs boson masses for Hi (red) and H 2 
(blue) as a function of ip K . Right: Tree-level (dashed) and one- loop (full) mass Mh 3 as a function of tp K . 
The exclusion region due to LEP, Tevatron and LHC data is shown as grey area, the region with the SM-like 
Higgs boson not being compatible with an excess of data around 125 GeV as dashed area. 
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Figure 2: The amount of CP violation r^p for Hi (i — 1,2,3) as a function of ip K (left). The amount of 
CP-even singlet component (Ti-a) 2 as a function of ip K (right). 



CP-conserving limit the one-loop masses of i?i 2 are Mh 1 = 119.4 GeV and Mh 2 = 120.7 GeV, 
with Hi being CP-even, cf. Fig. [2] (left), but CP-even singlet-like, cf. Fig. [2] (right), such that it 
hardly couples to SM particles and cannot be excluded by the experimental searches. The heavier 
Higgs H2 is dominantly CP-odd singlet-like (not plotted here) and is not excluded by the LEP, 
Tevatron and LHC searches due to both its singlet and its CP-odd nature leading to a vanishing 
coupling to weak vector bosons, cf. Fig. [3| With increasing CP-violating phase the eigenstates Hi 
and H2 interchange their roles both with respect to their CP nature and their amount of CP-even 
singlet component, with the cross-over taking place at (p K ~ it/ '64. 

In order to get an estimate of the theoretical uncertainty due to the unknown higher-order 
corrections the one-loop corrections to the Higgs boson masses have been calculated with the top 
and bottom pole quark masses, Mt = 173.2 GeV and Mj, = 4.88 GeV, and compared to the 
results for the one-loop corrected masses evaluated with the running DR top and bottom quark 
masses m t ^ at the scale Q = y jmQ z mt R . For our scenario they amount to mt = 153.4 GeV and 
mf, = 2.55 GeV. The result is shown in Fig. |4j Whereas the slope of the curve hardly changes, the 
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Figure 3: The Hi coupling to V (V — Z,W) bosons squared (i = 1,2,3) normalized to the SM coupling, 
\9vvHi\ 2 , as a function of ip K . 
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Figure 4: The one-loop corrected mass of the SM-like Higgs H3 evaluated with the top and bottom running 
DR masses (full) and with the corresponding pole masses (dashed). 

absolute values of the corrections change and are more important for a higher top quark mass. The 
theoretical uncertainty due to the different quark mass renormalization schemes can conservatively 
be estimated to ~ 10%. 



4.1.2 No tree-level CP violation 

We now keep the CP- violating phases <p K and ip\ non-zero and vary them by the same amount, such 



that according to Eq. (19) we have no tree- level CP-violating phase ip y . In the one- loop corrections 
ip K ,ip\ enter separately so that CP violation is induced radiatively. Figure [5] (left) shows the one- 
loop corrected masses of the three lightest Higgs states Hi,2,3, Fig. [5] (right) compares the tree-level 
and one- loop corrected mass of the SM-like H3, both as a function of (p K . The tree- level mass 
shows no dependence on <p K as expected. The one-loop mass Mjj 3 changes by only ~ 3 GeV for 
(p K varying from to tt, and the loop-corrected masses for H\ i show almost no dependence on the 
CP-violating phase. The reason is that the dependence on the phase is due to the corrections from 
the stop sector which are the dominant contributions to the one-loop masses. The values of the 
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Figure 5: One-loop corrected Higgs boson masses Mjj t (i = 1, 2, 3) as a function of <p K 
(dashed) and one- loop corrected (full) mass for H3 as a function of ip K = ip\ (right). 
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stop masses change with the CP-violating phase. As H3 has the largest h u component and hence 
couples more strongly to the up-type quarks it shows a stronger dependence on ip K than Hi and 
H2. For Mffj (Mh 2 ) the mass corrections are of about 15 (11) GeV. With mass values around 120 
GeV they could lead to additional signals at the LHC if they were SM-like. However, due to the 
CP-odd nature of H2, cf. Fig. [6] (left), it hardly couples to weak vector bosons. And the CP-even 
singlet character of H%, compare with Fig. [6] (right), reduces its couplings to SM particles. These 
particles would therefore have considerably reduced signals at the LHC. The whole region over 
which (p K = (f\ are varied is hence still allowed by the LHC searches. 
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Figure 6: The amount of CP violation r^p for Hi (i = 1, 2, 3) as a function of tp K — <px (left). The amount 
of CP-even singlet component (Ha) 2 as a function of ip K = ip\ (right). 



4.1.3 Radiatively induced CP violation through the stop sector 

For completeness we investigate the case where only (fA t ^ 0. CP violation is thus only induced 
through loop corrections stemming from the stop sector. The one-loop corrected masses of #1,2,3 
are shown in Fig. [7] (left), the tree-level and one-loop corrected mass of the SM-like H% are displayed 
separately in Fig. [7] (right), both as function of <pA t - The tree-level masses of H\^ are increased by 
about 10-15 GeV and their one- loop masses of ~ 119.5 and 121 GeV, respectively, hardly show any 
dependence on ipA t . The SM-like #3 one- loop mass shows a small dependence varying by ~ 2 GeV 
for (fA t E [0, 7r], increasing the tree-level mass by 4-7 GeV. The reason is that #3 has the largest 
h u component so that the dominant one-loop corrections stemming from the stop loops contribute 
more importantly to the radiative corrections of the Higgs mass matrix elements of H3 than of H\ 
and H2. The CP-even singlet nature of H\, the CP-even character of H\ and #3 as well as the 
CP-odd one of H2 are hardly affected by a change in ipA t and are therefore not displayed here. As 
may have been expected, in this scenario loop-induced CP violation affects the phenomenology of 
the Higgs bosons less. Note, that the scenario is not excluded by LHC searches over the whole 
displayed phase range. 

In Fig. [8] we investigate the theoretical error due to unknown higher-order corrections by varying 
the renormalization scale from 500 GeV to half and twice the scale. The variation of the renormal- 
ization scale also changes the values of the input parameters and the running DR top and bottom 
mass. For higher scales they become smaller and hence also the one- loop corrections to the masses 
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Figure 7 
(dashed) 



One- loop corrected Higgs boson masses Mb ( (i — 1,2,3) as a function of <pA t (left). Tree- level 
and one-loop (full) mass Mh 3 as a function of ipA t (right). 



decrease. The residual theoretical uncertainty can be estimated to about 4%. We also checked the 
theoretical uncertainty due to the different quark mass renormalization schemes and found them 
to be of ~ 10%, hence of the same order as in the scenario studied in Sect. 4.1 
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Figure 8: One-loop corrected Higgs boson masses Mh 3 as a function of tpA t for three different renormalization 
scales, fi ren = 250 (blue/long-dashed), 500 (green/short-dashed) and 1000 GeV (red/dotted). 



In summary, the discussion of the various scenarios has shown that the impact of the CP- 
violating phase is crucial for the validity of the model. While a certain parameter set can still 
accommodate the experimental results for vanishing CP violation it may be invalidated by non- 
vanishing CP phases. Turning this around, the experimental results will be useful to pin down 
the allowed amount of CP violation. The latter can arise from tree-level CP-violating phases in 
the Higgs sector or be radiatively induced. In the latter case the effects are found to be less 
pronounced. To get reliable predictions, the one-loop corrections have to be included as they 
not only considerably change the absolute mass values but also the singlet and CP-nature of the 
individual Higgs bosons as compared to the tree- level quantities. 
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4.2 Scenario with SM-like H\ or H 2 



The parameter set for this scenario, where, depending on the CP-violating phase, either Hi or H2 
is SM-like, is given by 

I A| = 0.65 , |k|=0.25, tan = 3, M H ±=619GeV, |A K | = 18GeV, \fi\ = 199 GeV 

\A b \ = 971 GeV , \A t \ = 1143 GeV , M 1 = 105 GeV , M 2 = 200 GeV , M 3 = 600 GeV . 

(116) 



The signs of the tree-level CP-violating phases Eqs. (|28|) and (|29|) are 

+1 , 



sign cos cp a 



sign cos ip z 



-1 . 



(117) 



The renormalization scale has been set to /J. ren — 650 GeV. The left- and right-handed soft SUSY 
breaking mass parameters in the stop sector m,Q 3 = 642 GeV and m tR = 632 GeV lead to mr = 
514 GeV and = 768 GeV. The low value of A respects the bounds imposed by unitarity 
We allow for tree- level CP violation by choosing ip K 7^ 0. The remaining complex phases are all set 



to zero, except for ipA x an d <PA K which follow from the tadpole conditions Eqs. (15), fll7h. 



The tree- level and and one-loop corrected masses of Hi and H2 are shown in Fig. [9] (left), 
as a function of cp K . Beyond ip K ~ 0.17T the tadpole conditions are not fulfilled any more. The 
corresponding couplings squared to weak vector bosons are plotted in Fig. [9] (right). The amount 
of CP violation of the three lightest Higgs bosons #1,2,3 and their CP-even singlet component are 



displayed in Fig. 10 (left) and (right), respectively, as a function of tp K . As can be inferred from 
the Figures, in the limit of the real NMSSM Hi is CP-even and has SM-like couplings while H2 is 
CP-odd. The heavier H3 is CP-even over the whole <p K range. The CP-even singlet components 



of Hi and H2 vanish, cf. Fig. 10 (right). However, H2 is CP-odd singlet-like. This is reflected in 
the couplings of Hi and H2 to the weak vector bosons. The one- loop corrections for <p K = shift 
the Hi mass from 99 to about 122 GeV so that its mass is compatible with the excess observed 
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Figure 9: Left: Tree-level (dashed) and one-loop corrected (full) Higgs boson masses as a function of tp K 
for H\ (red) and H 2 (blue). Right: The Hi (red) and H 2 (blue) Higgs couplings squared to two V bosons 
(V = W, Z) as a function of (p K at tree-level (dashed) and at one-loop (full). The exclusion region due to 
LEP, Tevatron and LHC data is shown as grey area, the region with the SM-like Higgs boson not being 
compatible with an excess of data around 125 GeV as dashed area. 
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at the LHC. The mass of the second Higgs boson H2 is increased by about 3 GeV to 126 GeV 
and could have been observed at the LHC, if its coupling to gauge bosons were not suppressed due 
to its CP-odd character so that it is not excluded by the present experimental constraints from 
the LHC. Furthermore, for ip K = the total significance concerning H\ is compatible with LHC 



searches according to our criteria Eq. (110). 



The scenario is interesting because with increasing cp K the CP character of H\p, changes rapidly 
(with a cross-over at (p K ~ 37r/64 where H\ t 2 interchange their roles with Hi being more CP-odd 
like and H2 more CP-even like). This dependence on the CP violating phase is at one- loop more 
pronounced than at tree-level and makes that already beyond (p K ~ 0.006 ir H\ cannot fulfill the 
role of the SM-like Higgs boson any more as its couplings deviate too much from the SM case 



to fulfill the requirement of Eq. (110). On the other hand, the if 2 couplings are not yet SM-like 
and once this is the case H2 is already too heavy to be compatible with LHC searches. Hence, 
above ip K ~ 0.006 ir the scenario does not comply with the criteria of Eq. (110) any more and it is 
excluded as indicated by the dashed region. Beyond (p K f» 0.052 tt the grey region shows that the 
searches at LEP invalidate this parameter choice due to the LEP limit on Hi in Zbb. Therefore a 
large portion of this scenario is likely to be excluded, constraining (p K to be almost zero and hence 
a real NMSSM. 

This scenario illustrates particularly well the importance of the one-loop corrections and the 
impact on the restriction of a possible CP-violating phase. Firstly, the one-loop corrections are 
crucial to shift the mass of the SM-like Higgs boson to a mass value which is compatible with 
the excess observed at the LHC. However, the one-loop corrections also amplify the dependence 
on the CP-violating phase of both the Higgs masses and in particular the mixing matrix elements 
and hence the coupling to the weak vector bosons. Neglecting for the moment for the sake of this 
discussion the fact that at tree-level Hi does not fulfill the mass constraint, the restriction of the 
CP-violating phase due to deviations from the SM significance would be less severe at tree-level 
than at one-loop level due to the smooth tree-level dependence on the CP-violating phase. The 
one-loop corrections are hence crucial to correctly define parameter scenarios which are compatible 
with present LHC searches and to derive the correct exclusion limits for scenarios dropping out of 
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Figure 10: The amount of CP violation r^p for Hi (i — 1,2,3) as a function of <p K (left). The amount of 
CP-even singlet component (7^3) 2 as a function of ip K (right). 
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this constraint. 

We close this subsection with the discussion of the CP-even H%, not shown explicitly in all plots. 
It is dominantly CP-even singlet-like. Its tree-level mass of ~ 148 — 152 GeV for {p K increasing 
from to ~ 37r/16 receives one-loop corrections of ~ 3 — 4 GeV. The one-loop corrections show the 
same dependence on tp K as the tree-level mass. Due to its singlet character at present it cannot be 
excluded by LHC searches. 

5 Summary and Conclusions 

We have calculated the one-loop corrections to the neutral Higgs bosons in the CP-violating 
NMSSM by applying a mixed renormalization scheme where part of the parameters are renormalized 
on-shell while tan/3, v 8 , A, k, A k and the CP-violating phases are renormalized in the DR scheme. 
We have in general allowed for tree- level CP violation due to non- vanishing phases (p u ,<p s ,(p K and 
(px, and for loop induced CP violation from the stop sector due to a non-zero phase <pA t - Several 
scenarios have been investigated which start from parameter sets that are compatible with the 
experimental Higgs searches in the limit of the real NMSSM, subsequently CP violation is turned 
on. As expected the dependence of the one-loop corrected Higgs masses and mixing matrix ele- 
ments on the CP- violating phase turned out to be more pronounced for tree-level CP violation than 
for radiatively induced CP violation. The loop corrections were found to considerably change the 
masses and mixing angles with crucial implications for the Higgs phenomenology at the LHC. As 
it is well known in the MSSM and the real NMSSM we also found that a scenario may be excluded 
at tree-level, whereas it is compatible with LHC searches at one-loop. Of special interest is the 
dependence on the CP-violating phase. It may be rather smooth at tree-level but more pronounced 
at one-loop so that at one-loop the CP-violating phase under investigation may be much more 
restricted than at tree-level due to possible non-compatibility with the experiments. Therefore, 
in order to correctly define viable scenarios and pin down allowed parameter ranges the one-loop 
corrections are indispensable. We also investigated the theoretical error due to the unknown higher 
order corrections by applying an on-shell and a DR renormalization scheme for the top and bottom 
quark mass and by varying the renormalization scale between half and twice its value. The theo- 
retical error of the one-loop corrected Higgs masses can be conservatively estimated to be about 
10%. 
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A Relations between Original and Physical Parameters 



For the transformation of the Lagrangian from the original parameters to the physical ones the 
following relations are used: 
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where n x and n z can be zero or one in case of two solutions of the tadpole condition, Eqs. (15) 



and (17), and zero if there exists only a single one. Here, sign x and sign 2 are the sign of the 



corresponding arcsine evaluated in the interval [— n, w), respectively. 
B Higgs Boson Mass Matrix 

In this section we list the Higgs boson mass matrix elements in a form needed as starting point for 
the renormalization procedure in the basis $ = (hd, h u , h s , A,a s ,G) T '. This basis is obtained by 
transforming the original basis (f> = (hd, h u , h s , a u , ad, a s ) T with the matrix 
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(126) 



The angle ft n is chosen such that the Goldstone boson field (with zero mass eigenvalue) is extracted 
and, at tree-level, coincides with the angle ft defined via the ratio of the vacuum expectation values, 

ftn = ft- 
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The mass matrix elements of the CP-even part Mhh, cf. Eq. (21), are 

M h d h d = - M$y]s} + M 2 c 2 + 9 J C f% [(1 + 2tptp B )t hd - t p t hu ] 

+ 2|A| 2 M^%4, (127) 



M hd hu = ~ [-^ - Myy + Mf] w + — fe [t hu + tpt 2 0B t hd ] 

„2 



-!- |A| 2 M^^f ^ , (128) 



M huhu = [^M " M^]c 2 + Mh} + [(2 + tpt PB )t hu - tp B t hd 



+ 2\A\ 2 M^4, (129) 

^hdft. - - [72 M W\ — + —3 [thu + t^pE^A 

c A/3 ttUs u s c A/3 

+ |A|il%^ s [2|A| C/3 - \K\s^ y ] - 4|A|2M ^^ g ' C/3 , (130) 
e e v s 

r2 ^„..„„ __„__ AA 



M h ,„ = - [ V - Mi=] M ^f^ + ^K + .AU 

+ \\\M w S ^v s [2\\\sp-\ K \cpc ipy ] - WJwAA (131) 

e e J f s 

M h , ft , = [-2 M w \ -2-2 —^2 ^hu + Wp B th 4 \ + — 



c A/3 e V s ev s c A0 

+ |A|M^%^ [2|A|M^%s 2(S - |k|^%J + 2\k\ 2 v 2 s + -U^H^c^ , (132) 

where A/3 = (3 — /3b and (3b = fic = fin- The mixing angle of the charged Higgs bosons, /3 C , and 
the mixing angle (3 n , needed for extracting the Goldstone boson, coincide and no discrimination 
between these two mixing angles is done in the formulae. The angles <p x , <p y and (p z have been 
defined in Eqs. Jl8)-(|20). 

The mass matrix elements corresponding to the CP-odd components of the Higgs boson mass 
matrix are given as 

s 2 

M AA = M 2 H± - M^c 2 Ap + 2\\\ 2 M^-^-c? Ap , (133) 
M Aas = [M 2 H± - M 2 w cU - ^ [t K + t,t 2 B t hd ] 



\X\M w S ^^[2\X\M 2 v S ^s 2 p-3\K\v 2 s c ipy ] , (134) 



30 



M as 



e v s c A/3 



ev s c A/3 



+ \X\M^ S -^^[2\X\M^ S -^s 2( s + 3\K\v 2 s c Vy ] - ^=\A k \\k\v s c Vz , 

6 V s 6 v ^ 



Mag = [M 2 H± - M^cl^tAp + 



ecp E 



2MiyS6» w CA/3 



(135) 
(136) 



et> s c 



A3 



A/3 



+ |A|A% 



S0 w SA/3 



ev. 



[2\\\M' 2 v \s 2f) - MvW y ] , 



M GG = [M|± - Af&c^,]*^ + 



eC/3-2/3 B 



(137) 



[t hd - tp- 2 p B t hu ] + 2\\\ 2 M%r-^slp . (138) 



Finally, the mass matrix elements describing the mixing between the CP-even and the CP-odd 
components can be expressed as 



M ha U gT 



( 2MZLs^ d ±ta d +3\K\\\\M W S -^V sS pS Vy - 2M Zl: w s, ta d \ 



es 

V M hsA 



2M w sa„ r s« ta d vlt p ta d + 3 M I M M W °Y v s c f3S l p y 2 M w sg w sp ta d 

M hsas 



ec P B 

M hsG ) 



(139) 



with 



M haA = ^-t ad ~ \k\\\\M W ^Wv s CA(3S^ , 

v s s/s e 

Jl-f 2 + ^ M wse w rCp . , ma , s e -, 
Mft s0s = — t as -V^aa + |k| \X\M W —^-8 2 pS Vy \ , 



v s s/3 e 
C Chargino and Neutralino Self-Energies 



(140) 

(141) 
(142) 



In this section the expressions for the renormalized self-energies of the charginos and neutralinos 
are listed. The different parts of the renormalized chargino self-energy decomposed according to 
Eq. (73) can be written as (i,j=l,2) 



(143) 
(144) 





(P f, )]« = 




(P 2 ) 


]« + 


2^*(^ + ^)C/ T ].. 5 




(^)],- = 




(P 2 ) 


]« + 


^[n^f +szT)v% , 




(^)],= 




(P 2 ) 


Lj 


^m x ±([C/*5Zg?7 T ] ifc (5 fej + ^[y^f^fe,-) - 




(P%i = 


PR 


(P 2 ) 




\m x ±([V5Z^V% 5 kj + 5 ifc [[T<5Z£*[/ T ] fej ) 
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where, for the renormalization procedure, the chargino spinors as given in Eq. ^ and the 2x2 
chargino mass matrix Mc are replaced by 

^+^{l+Uz c L y + L , (147) 
1>R -+ (*+ \$ Z r)^R > with SZg = (y Z ^ ) and X = L,R , (148) 

and 



M c ^t M c + SM C , (149) 



respectively, where 8Mq is given in Eq. (78). 

The various parts of the decomposed renormalized neutralino self-energy can be expressed as 
(i, j=l,...,5) 

[^(p 2 )]* = + \ W*V zS + sz^V^ia , (iso) 

K° ^ a = K° &)] a + \ WW" + «^ ^ « > ( 151 ) 

- [N*5M N N% , (152) 

[£JV)L, = KoCp 2 )]^ - ^m x o([A/'52' iV *j'v rT ]ifc + S i k[M5Z N "'M T ]kj) 

- [M5M ] N M T ] i3 . (153) 

For the renormalization procedure, the neutralino spinor defined in Eq. Q has been replaced by 
i/>° ^ (l+l5Z N )l>° with 6Z N = diag(6Z N ,6Z?,6Zg,8Z?,6Zg) (154) 



M N -4 M w + 5Mjy ■ (155) 



2 

and the neutralino mass matrix by 



The matrix elements of JM^r can be found in Eqs. (86)-(96) 
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